Mon. Not. R. Astron. Soc. 000, 000-000 (0000) Printed 2 April 2013 (MN Mg^ style file v2.2) 

The metallicity - redshift relations for emission-line SDSS galaxies: 
examination of the dependence on the star formation rate 



en 

o 



L.S. Pilyugini'2, M.A. Lara-Lopez^ E.K. Grebel^, C. Kehrig'*, LA. Zinchenko^ 
A.R. Lopez-Sanchez^'^, J.M. Vilchez^, L. Mattsson^ 

' Main Astronomical Observatory of National Academy of Sciences of Ukraine, 27 Zabolotnogo str, 03680 Kiev, Ukraine 

^ Astronomisches Rechen-Institut, Zentrum fur Astronomic der Universitdt Heidelberg, Monchhofstr 12-14, 69120 Heidelberg, Germany 

Australian Astronomical Observatory, PO BOX 915, North Ryde, NSW 1670, Australia 
* Instituto de Astrofisica de Andalucia, CSIC, Apdo, 3004, 18080 Granada, Spain 

Department of Physics and Astronomy, Macquarie University, NSW 2109, Australia 
^ DARK Cosmology Centre, Niels Bohr Institute, University of Copenhagen, Juliane Maries Vej 30, DK-2100, Copenhagen 0, Denmark 



cn 



Accepted 2013 March 26. Received 2013 March 03; in original form 2013 January 1 1 



o 
u 

o 



> 

a^ 

o 

^' 

o 
m 



^ 



ABSTRACT 

We analyse the oxygen abundance and specific star formation rates (sSFR) variations with 
redshift in star-forming SDSS galaxies of different masses. We find that the maximum value 
of the sSFR, sS FR„a^-, decreases when the stellar mass, Ms, of a galaxy increases, and de- 
creases with decreasing of redshift. The sSFR^ax can exceed the time-averaged sSFR by 
about an order of magnitude for massive galaxies. The metallicity - redshift relations for sub- 
samples of galaxies with sSFR - sSFR„,ax and with sSFR - Q.lxsSFR„ax coincide for mas- 
sive (log(Ms/Mo) > 10.5, with stellar mass Ms in solar units) galaxies and differ for low-mass 
galaxies. This suggests that there is no correlation between oxygen abundance and sSFR in 
massive galaxies and that the oxygen abundance correlates with the sS FR in low-mass galax- 
ies. We find evidence in favour of that the irregular galaxies show, on average, higher sSFR 
and lower oxygen abundances than the spiral galaxies of similar masses and that the mass - 
metallicity relation for spiral galaxies differs slightly from that for irregular galaxies. The fact 
that our sample of low-mass galaxies is the mixture of spiral and irregular galaxies can be 
responsible for the dependence of the metallicity - redshift relation on the sS FR observed for 
the low-mass SDSS galaxies. The mass - metallicity and luminosity - metallicity relations 
obtained for irregular SDSS galaxies agree with corresponding relations for nearby irregular 
galaxies with direct abundance determinations. We find that the aperture effect does not make 
a significant contribution to the redshift variation of oxygen abundances in SDSS galaxies. 
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1 INTRODUCTION 

It is well established that the oxygen abundance correlates 
with galaxy mass (or luminosity), O/H = f(Ms), in the sense 
that the higher the ga laxy mass (luminosity), the hi gher the 
heavy element content JLequeux et al. 1979; Gamett & Shields I 



ivy 

si. 



19871; ISkillman et aT] 



1989; IVila-Costas & Edmunds I 



Zaritskv et al. I Il994l; iPilyuginI l2001bl ; iGuseva et al 



1992; 



20091 ; 



Lopez-Sanchez & Estebanll201CI , among others). In recent years, 



the number of available spectra of emission-line galaxies has 
increased dramatically due to the completion of several large 
spectral surveys, in part icular the Sloan Digital Sky Survey 
(SDSS) JYork et al.l 120000 . Measurements of emission lines in 
SDSS spectra have been carried out for abundance determina- 
tions and investigation of the luminosity - metallicity relation 
(e.g. iTremonti et al. 1 120041) . It is discussed whether the mass 



- metallicity relation can be caused by the galaxy-downsizing 
effect, where star formation ceases in high-mass galaxies at 
earlier times and shifts to lower-mass galaxies at later epochs 
JCowie et alj 119961) . The study of oxygen abundance evolution 
with redshift in emission-line SDSS galaxies pro vides strong 
evidence in favour of the galaxy-downsizing effect |Thuan et al. I 
I2OIOI ; iLara-Lopez et al. l2010al ; IPilvugin & Thuan II2OI ll) . 

It is important to note there are different absolute O/H 
scales. The empirical metallicity scale is defined by the direct 
meth od (Tg method) a nd the empirica l calibrations based on it 
(e.g. JPilvuginl 12001 j; IPettini & PageT] [2004; Pilvugin & ThuanI 
1 20051 ; IPilvugin et al. I uOlOh . Metallicities deri ved using calibra- 



tions based on photoionisatio n models (e.g. 



iKewlev & Dopit a 2002; Kobul nickv & Kewlev 



McGaugh|[T99il ; 
2004 ) tend to be. 



systematically higher (~0.4 dex or ev en more) th an those derived 
using the direct method (see reviews bv lLopez-Sanchez & EstebanI 
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uOlu iLopez-Sanchez et al. Il2012h . Therefore, there are large dis- 
crepancies between oxygen abundances in the SDSS galaxies ob- 
tained in different works using different calibrations. 

It has been argued that the mass - metallicity relation is due to 
a more general relation in the 3-dimensional space formed by the 
stellar mass, gas-phase oxygen abundance and the star formation 
rate. This 3-d imensional structure is ref erred to as a Fundamental 
Plane (FP) bv iLara-Lopez et al. I ( l2010bl) , suggesting that the stellar 
mass can be expressed as a functi on of the S FR and 0/H , i.e. Mj = 
f(SFR,0/ll). On the other hand. lMannucci et a l. ' flOlif) describes 
this 3-dimensional space as a surface, the Fundamental Metallicity 
Relation (FMR) in which the O/H = f(M„SFR). For a further dis- 
cussio n of the 3-dimensional stru cture of this space see lYates et al. I 
( I20 1 2h : iL^a-Lopez et alTI ( |20 1 3h . 

It is also known that the mass - metallicity relation changes 
with redshift z- Thus, one can expect that in general the oxygen 
abundance in a galaxy coiTelates with three parameters, i.e., O/H 
= /(A/5, S FR, z)- Study of the redshift evolution of oxygen abun- 
dances in galaxies of different masses may tell us something about 
the evolution of galaxies and provide important constrains on the 
models for the chemical evolution of galaxies. In this paper, we 
will examine whether the metallicity - redshift relations for SDSS 
galaxies of different masses are dependent on the star formation 
rate. 

The paper is structured as follows. Selection criteria for the 
sample of the SDSS galaxies with reliable abundance determina- 
tions are reported in Section 2. In Section 3 we establish the redshift 
evolution of the maximum star formation rates in galaxies of differ- 
ent masses as well the redshift evolutions of oxygen abundances in 
galaxies of different masses and sS FR. The properties of iiTegular 
and low-mass spiral galaxies are compared in Section 4. In section 
5 we discuss whether the aperture effect makes a significant contri- 
bution to the redshift variation of oxygen abundances in the SDSS 
galaxies. We also estimate the time-averaged values of the sS FR in 
massive galaxies. The conclusions are given in Section 6. 

Throughout the paper, we will be using the following nota- 
tions for the line fluxes, 



are described in lBrinchmann et al. I (120041) : iTremonti et al. I ( 120040 
and various other publications by the same authors. 

In a first step, we extract all emission-line objects with 
measured fluxes in the Hj8, Ha, [Oii]/t/l3727,3729, [O iii].15007, 
[Nn]i6584, [Sii]i6717 and [Sii]/16731 emission lines. The hydro- 
gen, oxygen, nitrogen and sulfur lines will serve to estimate oxy- 
gen and nitrogen abundances relative to hydrogen. The wavelength 
range of the SDSS spectra is 3800 - 9300 A so that for nearby 
galaxies with redshift z < 0.023, the [On]i3727+/13729 emission 
line is out of that range. Therefore all SDSS galaxies with z < 0.023 
were excluded. In the other end, distant galaxies with redshift z 
> 0.33, the [Sii]/16717-t-/16731 emission line is out of that range. 
In general, reliable estimates of the oxygen abundance in a SDSS 
galaxy can be fo und even if the [O ii]/l3727-l-/l3729 emission line 
is no t available jPilvugin & Thuan 1 120071 : IPilvugin & Mattsson I 
l201lh or if the [Sii]/ t6717-l-/l6731 emission line is not measured 
I Pilyugin et al. I2OI0I) . However, our strategy to construct a sample 
of the SDSS galaxies with rather reliable oxygen abundances re- 
quires both these lines. Thus, all galaxies in our total sample have 
redshifts greater than ~ 0.023 and lower than ~ 0.33. The redshift 
z and stellar mass M^ of each galaxy were also taken from the 
MPA/JHU catalogs. 

The emission-line fluxes are then corrected for int erstellar red- 
dening using the theoretical llafH ratio and the Whi tford I ( Il958h 
interstellar reddening law due to llzotov et al. I ( I1994J) . We adopt 
case B recombination with a density of 100 cm"^, and an elec- 
tronic temper ature of 10"^ K, which yields a theoretical ratio of 
Hq'/Hj8=2.86 tosterbrock & Ferland II2OO6I) . In several cases, the 
derived value of the extinction C(H/3) is negative and has then been 
set to zero. 

The [Oiii]/i4959 line is required to obtain the Rj, value. 
However, this line is measured with much less accuracy than 
[Oiii]i5007 line. Since the [Oiii]/I5007 and A4959 lines originate 
from transitions from the same energy level, their flux ratio is due 
only to the transitio n probability ratio, which is very close to 3 
dStorev & Zeippen 112000.) . Therefore, the value of R3 can be es- 
timated as 



^2 = -'[OII],13727+,13729/-'h/3, 

^^2 = -'[Nll]/i6548+,i6584/-'n8, 

'5'2 = -'[SIIli6717+i673l/-'n8> 

R3 = -'[Om],M959+,15007/-^H/3- 



(1) 

(2) 
(3) 
(4) 



R3 = 1.33x[Oni].l5007/H/? 



(5) 



The stellar masses of galaxies Ms are expressed in the solar mass 
units. 



2 DATA BASE 

2.1 A sample of SDSS galaxies 

Our study is based on SDSS data. Line flux measurements in SDSS 
spectra have been caiTied out by several groups. We use here the 
data made publicly available by the MPA/JHU groujiJ. These data 
catalogs give line flux measurements, redshifts and various other 
derived physical properties such as stellar masses for a large sample 
of SDSS galaxies. The techniques used to construct the catalogues 



' The catalogs ai'e available in the SDSS DR9 database (tables galSpecEx- 
tra, galSpecInfo and galSpecLine) at^http://www.sdss3.org/dr9/, 



Similarly, the value of A'2 is estimated without the line [Nn]/i6548. 
The [Nn]/i6584 and A6548 lines also originate from transitions 
from the same energy level a nd the transition probabil ity ratio for 
those lines is again close to 3 dStorev & Zeippen 11200(1) . The value 
of Ni can therefore be estimated as 



N2= 1.33x[Nii].16584/II;e 



(6) 



The spectroscopic data so assembled form the basis of the present 
study. 

The intensities of strong, easily measured lines can be used 
to separate different types of em ission-line objects ac cording to 
their main excitation mechanism. [Baldwin et al. I ( 1198 ih proposed 
a diagram (the so called BPT classification diagram) where the 
excitation properties of Hii regions are characterised by plot- 
ting the low-excitation [Nii]/i6584/HQ' line ratio against the high- 
excitation [O in]/l5007/II/J line ratio. The exact location of the di- 
viding line between Hn regions and AGNs is still controversial , 
however (see, e.g., Kewlev et al. I uOOll : iKauffmann et al. I l2003l : 
St asihska et al. 2006 i ). Her e we adopt the dividing line suggested 
bv lKauffmann et al. I ( I2OO3I) 



log([Oin]i5007/II/?): 



0.61 



log([Nii]^6584/HQ') - 0.05 



+ 1.3, 



(7) 
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Figure 1. Comparison of differences between oxygen (nitrogen) abundances determined in different ways for the sample of reference H n regions (E sample). 



which separates objects with H ii spectra from those containing an 
AGN. By this formula, the suspected AGNs have been excluded 
from further consideration. 



2.2 Abundances 

A new method (labelled "the C method") for oxygen and ni- 
trogen abundance de termination has recently been suggested 
dPilvugin et al. luOlZ c?). It is based on the standard assumption that 
Hn regions with similar strong-line intensities have similar phys- 
ical properties and abundances. From a sample of reference H ii 
regions with well-measured (Je-based) abundances we choose a 
counterpart for the considered H ii region by comparison of four 
combinations of strong-line intensities: P = Ri/iRi + Ri) (exci- 
tation parameter), logR^, \og(N2/R2), and \og{S2/R2) The oxygen 
and nitrogen abundances, as well as the electron temperature in the 
studied Hn region may then be assumed to be the same as in its 
counterpart. To get more reliable abundances, we select a number 
of reference H ii regions with abundances near those in the coun- 
terpart H n region and then derive the abundance in the studied H n 
region through extra-/interpolation. 

It has been suggested that reliable oxygen and nitro- 
gen abundances in Hn region can be derive d even if the 
[S a] A67 n +A6731 emission line is not measured (Pilvu ginet al. I 
UOIOi) or if the [O n]/l3727-l-/l3729 emission line is not available 



jPilvugin & Mattsson luOl ih . Hence, the C method can be modi- 
fied in the following way. To find the counterpart for the H n region 
under study, we will compare not the four values that are expressed 
in terms of the strong line intensities but two sets of three values: 
I) logi?3, P and log(N2/R2) and 2) \0gR3, logA'a and log(A'2/52). 
The oxygen and nitrogen abundances determined with the help of 
the counterparts selected with the former set of values will be re- 
ferred to as (0/H)coN and (N/H)cof,5 while the oxygen and nitro- 
gen abundances determined with the help of the counterparts se- 
lected with the latter set of values will be referred to as (0/H)cns 
and ( N/H)cm„ . Using the co llected data with T^-based abundances 
from Pilvugin et al. I ( l2012cl) we select a sample of reference H 11 re- 
gions for which all the absolute differences for oxygen abundances 
(0/H)roN ~ (0/H)r, and (0/H)c^,s - (0/H)t„ and for nitrogen abun- 



dances (N/H),;- 



(N/H)r, and (N/H), 



Cns 



(N/H)r, are less O.I 



dex. This sample of reference Hii regions contains 233 objects 
and will be used for abundance determinations in the following. 
This sample will be referred to as E sample (etalon sample) below. 
Hence, here we use the empirical metallicity scale defined by the 
H II regions with abundances derived through the direct method (T^ 
method). Comparison between the oxygen abundances (O/H)^^,^ 
and (0/H)cq„ and nitrogen abundances (N/H)c„j and (N/H)fQ„) 
provides a selection criterion for star- forming galaxies with accu- 
rate line fluxes measurements (see below). 

The SDSS spectra of distant galaxies are closer to global spec- 
tra of galaxies than spectra of single H n regions. IPilvugin et al. I 
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( l2012bf) have Monte Carlo simulated global spectra of composite 
nebulae as a mix of spectra of individual components, based on 
spectra of well-studied H n regions in nearby galaxies. Abundance 
analysis of the artificial composite nebulae yielded the conclusion 
that oxygen and nitrogen abundances determined using the ON and 
NS calibrations are in good agreement with each other and near 
(within ~0.2 dex) the mean Hj8 luminosity- weighted value of oxy- 
gen and nitrogen abundances of individual components of the com- 
posite nebula. 

2.3 Selection criteria 

Our study is based on SDSS spectra. To study the redshift evolu- 
tion of oxygen and nitrogen abundances, accurate oxygen and nitro- 
gen abundance determinations are of course a prerequisite. Abun- 
dances in H II regions determined from SDSS spectra suffer from 
the following problem, however. Line fluxes in SDSS spectra are 
measured by an automatic procedure, which inevitably introduces 
large flux errors for some objects. We therefore wish to exclude 
them from further consideration. A method to recognise such ob- 
jects relies on the fact that there exist a fundamental Hii region 
sequence. The emission line properties of photoionised H ii regions 
are governed by their heavy element content and by the electron 
temperature distribution within each photoionised nebula. The lat- 
ter is controlled by the ionising star cluster's spectral energy dis- 
tribution and by the chemical composition of the H ii region. The 
evolution of a giant extragalactic H ii region associated with a star 
cluster is thus caused by a gradual change over time (due to stellar 
evolution) of the integrated stellar energy distribution. This phe- 
nomenon has been studied in numerous investigations JStasiiiska 



197 lll98d:lMcCall et al. ll98l:lDopi ta & Eva ns I1986| ; IMov et al. 



2001 



Stasiiiska & Izotov I l2003l : iDopita et al. l200d . among oth 



ers). The general conclusion from those studies is that Hii regions 
ionised by star clusters form a well-defined fundamental sequence 
in different emission-line diagnostic diagrams. The existence of 
such a fundamental sequence forms the basis for studying vari- 
ous p roperties of extragalactic Hii regions. First. Baldwin et al. I 
( 1198 If) have suggested the position of an object in some well- 
chosen emission-line diagrams can be used to separate Hii re- 
gions ionised b y star clusters from oth er types of emissio n-line 
objects. Second. [Pagel et al. I ( 1197^ and lAlloin et al. I ( Il979l) have 
suggested that the positions of H ii regions in some emission-line 
diagrams can b e calibrated in terms of their ox ygen abundance s. 

Following lThuan et al. I <201Cl) and .Mannucci et al. I (2010|) we 
adopt here a simple method to recognise star-forming galaxies with 
accurate line fluxes measurements. It is based on the idea that if i) 
an object belongs to the fundamental H ii region sequence, and ii) 
its line fluxes are measured accurately, then the different relations 
between the line fluxes and the physical characteristic of H ii re- 
gions, based on different emission lines, should yield similar phys- 
ical characteristics (such as electron temperatures and abundances) 
of that object. 

Fig. [T] show comparisons of differences between values of 
oxygen (nitrogen) abundances determined in various ways for the 
E sample of reference Hii regions. Fig.[T]also shows that the cor- 
relation between differences in oxygen abundances determined in 
different ways is weak. A correlation between differences in nitro- 
gen abundances determined in different ways is more evident. Al- 
though these correlations are not very tight one can expect that the 
conditions |log(0/H)c„, -log(0/H)coJ < 0.05 and |log(N/H)c„, - 
log(N/H)co„| < 0.10 allow to select the SDSS objects with more 
or less reliable abundances. Thus, to select star-forming galaxies 



with accurate line fluxes measurements, we require that, for each 
galaxy, the oxygen (0/H)c„s and (0/H)co„ and nitrogen (N/H)c„5 
and (N/H)co^ abundances agree. 

It should be noted that the C method was developed for H n re- 
gions in the low-density regime. Therefore we include only objects 
with reliable electron densities n^. The electron density is estimated 
from the sulfur line ratio [S ii]/16717/[S ii]/l6731. According to the 
five-level atom solution for the S+ ion, the [Sii]/l6717/[Sii]/l6731 
is a reliable indicator of the electron density in the interval from n,. 



~ 10^ cm-3 to n,~W cm"^ 

log n, = 8.448 - 15.101x -I- 14.419^^ - 5.115x3 



(8) 



where x = [Sii]/16717/[S ii]/16731. We exclude objects with 
[Sii].i6717/[Sii]/16731 < 1.2, i.e., those with n, > 170 cm-\ At 
densities n, < 10^ cm"^ the ratio [Sii]/16717/[Sii]/16731 is al- 
most independent on the electron density and reaches the limit 
value 1.43. The measured ratio [S n]/i6717/[S n]/l6731 in some ob- 
jects is in excess of this limit value, that can be considered as 
evidence in favour of that the line measurements in the SDSS 
spectra of this object are not accurate. Therefore the objects with 
[S ii]/i6717/[S ii]/16731 > 1.55 (expected error in the line measure- 
ments is larger ~10%) were also excluded. 

Another subset of galaxies were excluded from consideration 
using the following conditions. In some cases, H ii regions in giant 
spiral galaxies seem to be catalogued as dwarf galaxies of small 
dimensions. Thus, if a "galaxy" has a Petrosian radius less than 3 
kpc, this galaxy is then excluded from our list. Further, the aper- 
ture correction factor (see below) for some galaxies is < 1, i.e., the 
galaxy light within the fiber exceeds the total light of the galaxy. 
The galaxy light fraction within the fiber is less than 0.7 for the 
bulk of galaxies. Therefore, the objects with the galaxy light frac- 
tion within the fiber higher 0.8 were also excluded from further 
consideration. 

Our final sample contains 98.986 SDSS objects. The O/H - 
N/0 diagram can tell us something about the credibility of our se- 
lection criteria and about the reliability of the abundance determi- 
nations. Fig. [2] displays the N/O ~ O/H diagram. The gray (blue 
points) in the panel a) show the (X/H)q^^ abundances for the se- 
lected sample of the SDSS galaxies, the gray (blue points) in the 
panel b) show the (X/H)c^g abundances. The dark (black) points 
(in both panels) show the T^-based abundances for the sample of 
reference Hii regions (E sample). A quick look at Fig. |2] immedi- 
ately tells us that the selected SDSS galaxies, with very few ex- 
ceptions, are located within the area outlined by the the reference 
Hi: regions with T^-based abundances (as in the case (X/H)co^ and 
as in the case (X/H)c^j abundances in the SDSS galaxies). This 
suggests that the strategy outlined above allows us to construct a 
sample of SDSS galaxies with rather reliable oxygen and nitrogen 
abundances. 

Panels c) and d) in Fig. [2] show the N/O - O/H diagrams for 
the rejected SDSS objects. The positions of rejected objects oc- 
cupy a much larger area in the (N/0)co^ - (0/ii)coN (and (N/0)c„5 
- (0/H)c„j ) diagram than that outlined by the reference H ii re- 
gions, although some rejected objects are located within that area. 
The weak point of our selection criterion is that we lose some 
fraction of objects with reliable (0/H)q^^ abundances if only sul- 
fur [Sii]/16717 and [Sii]/16731 lines measurement involve eiTors, 
and some fraction of objects with reliable (0/H)c^,5 if only oxygen 
[O ii]/l/13727,3729 emission lines measurement involve errors. 

The resultant reduced sample is the basis for our study. For 
clarity and simplicity, only (0/H)co;^ abundances will be consid- 
ered in the following. 
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Figure 2. The N/0 - 0/H diagram. The gray (blue) points in panel a) show the (X/H)f ^^ abundances for the selected sample of the SDSS galaxies, the gray 
(blue) points in panel b) show the (X/H)c;^,s abundances. The gray (blue) points in panels c) and d) show the same for the rejected SDSS objects. The dark 
(black) points in the all panels show the r^-based abundances for the sample of reference H ii regions {E sample). (A colour version of this figure is available 
in the online version.) 



3 THE REDSHIFT VARIATIONS OF SFR AND OXYGEN 
ABUNDANCES IN SDSS GALAXIES 

3.1 Specific star formation rate 

The specific star formation rate sS FR is defined as the S FR per 
unit stellar mass in a galaxy, i.e., sS FR = S FR/Ms- In estimating 
S FRs from the Ha luminosity we adopt the calibration suggested 
bv lKennicuttl J1998h 



SFR (Moyr-'): 



(W) 



1.27x103*' 



where 



I^Ha ■ 



4nd'FZ, 



(9) 



(10) 



is the Ha luminosity. In equation dlOt d is the distance to a galaxy, 
and F™^ is the Ha flux corrected for extinction and aperture effect. 
The reddening-corrected HyS flux F^ is obtain ed from the ob- 
served F?^ flux using the relation JBlagrave et al. Il2007h 



log(F*/F^) = -Cus, 



(11) 



where C(IIy0) is the extinction coefficient. The generic reddening- 
corrected flux ratio is assumed to be F^JF^ = 2.86. 



In addition to the extinction correction, the Ha luminosity also 
requires an aperture correction to account for the fact that only a 
limited amount of emission from a galaxy is detected through the 3" 
diameter fiber We estimate the aperture correction factor A using 
the aperture correction method from ,Hopkins et al., ( 2003.) 



\Q-OA{rp„ra-rfii,e,) 



(12) 



where rfii,„ and r/j^,,,-,, are r-band fiber and Petrosian magnitudes. 
This correction is based on the assumption that the emission line 
flux can be traced overall by the r-band emission. The aperture- 



corrected Ha flux is F^°^ 



■-AxF2, ■ 



Distances to galaxies were calculated from 

cz 



(13) 



where d is the distance in Mpc, c the speed of light in km s ', and 
z the redshift. Hq is the Hubble constant, here as sumed to be equal 
to 72 (±8) km s"' Mpc ' (Freedman et al. 11200 ih . 
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Figure 3. Left panels. Redshift dependence of the sSFR in galaxies of different masses. Light-gray (blue) points ai'e individual galaxies. The solid line is the 
"sS FR„ax - redshift" relation, while the dashed line show the "0. 1 xsS FR„,ax - redshift" relation. Right panels. Redshift dependence of the oxygen abundance 
in galaxies of different masses and different sSFR. Dark-gray (red) points are galaxies with \og{sS FR)=\og(sS FR„,ax)±^).2. The solid line is the best fit to 
those data. The light-gray (blue) pints are galaxies with log(.vSFR)=log(0. lXiSFi?„ioi)±0.2, while the dashed line is the best fit to the data. (A colour version 
of this figure is available in the online version.) 
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3.2 The redshift variations of sS FRs and oxygen abundances 
in galaxies of different masses 

Studies of the redshift evolution of the 5 FRs in the SDSS galaxies 
encounter selection-effect problems in the sense that the weak hy- 
drogen emission line luminosity cannot be measured in the distant 
objects. The lower limit of the measured Hy8 luminosity at the red- 
shift 2=0.25 is three orders of magnitude higher tha n that at z=0.01 
(see, for example, Fig.7 in lPilvugin et al. I ( l2012ah ). To avoid this 
selection effect we will consider primarily the redshift evolution of 
the maximum value of the specific star formation rate sS FR„,^^ in 
galaxies of different masses. 

We consider six stellar-mass bins centered on the values logM, 
= 9.0, 9.5, 10.0, 10.5, 11.0, 11.5 with interval width +/- 0.2 dex. 
The left panels in Fig. [3] show the redshift variations of sS FRs in 
galaxies of different masses. Light-gray (blue) points are individ- 
ual galaxies. The solid line is the linear "sSFR„c,x - redshift" re- 
lation. This relation has been obtained as follows. As a first step, 
the average value of the sS FR has to be found. In a second step, 
the galaxies with sS FRs lower than the average value have to be 
rejected and a new average sS FR value has to be obtained. In a 
third step, the galaxies with the sS FRs lower than the new average 
value is rejected and a "sS FR - redshift" relation is then obtained. 
The fourth step is to reject the galaxies with the sS FRs below this 
relation and thus a final "sS FR„,ai - redshift" relation has been ob- 
tained. In the deriving the final "sS FR„ai - redshift" relation the 
objects with deviations exceeding the mean value of deviations by 
a factor three were excluded. Hence, the final "sS FR„,a^- - redshift" 
is determined using around 5-10% of galaxies in a fixed mass in- 
terval. In the first and second steps we use the average sS FR value 
rather than "sS FR - redshift" relation in order to exclude the ob- 
jects with low i'S FRs in an attempt to overcome the selection effect. 
The dashed lines in the left panels in Fig.|3]show the "0. IXiS FR^ax 
- redshift" relations. 

The left panels in Fig. [3] shows that the maximum specific 
star formation rate correlates with the stellar mass of a galaxy: the 
value of the sS FR,„ax decreases with increasing of the Ms . This 
nicely confirms the galaxy downsizing effect, where star formation 
ceases in high-mass galaxies at earlier times an d shifts to lower- 
mass galaxies at later epochs JCowie et al.l 19960 . 

The right panels in Fig. [3] show the redshift variations of the 
oxygen abundance in galaxies of different masses and different 
sS FRs. Dark-gray (red) points are galaxies with sS FRs close to 
sS FR„,„x, i.e., galaxies which are located within +/-0.2 dex along 
the solid line in the corresponding left panel in Fig. |3] The solid 
line is the linear best fit to those data. The light-gray (blue) pints are 
galaxies with sS FRs close to O.lxsS FR,„i,x, i.e., galaxies which are 
located within the band of width +/-0.2 dex along the dashed line 
in the corresponding left panel in Fig.|3] The dashed line is the best 
fit to those data. 

A closer look at Fig. [3] reveals that the metallicity - redshift 
relations for massive galaxies (logM, > 10.5) do not depend on the 
sS FR. While the sS FR changes by an order of magnitude (from 
sS FRi„„x to O.lxsS FR„„x), the metallicity - redshift relations for 
those galaxies are similar. This, we interpret as evidence that there 
is no correlation between oxygen abundance and sS FR in massive 
galaxies. The metallicity - redshift relation for low-mass galaxies 
(logM, < 10.5) on the other hand, depends on the sS FR: the galax- 
ies with highest sS FR show lowest oxygen abundances. 



3.3 Comparison to a previous analysis of SFR-dependence of 
the O/H - M, relation 

The mass - metallicity relation of galaxies as a function of (spe- 
cific) star formation rate has previously been considered in sev- 
eral different stu d ies JEUison et al.l2008l:lLara-L6Dez et al. luOlObl: 



iMannucci et al. 20 id: 



2013; Lara -Lopez et al. I 



Yates et al. 



2OI2I : [Andrews & Martini I 



20131) . It is difficult to directly compare 



those results with ours for two reasons. First, our method for abun- 
dance determination differs from that used in other studies. Second, 
a single mass - metallicity relation for galaxies from a rather large 
redshift range were constructed in the papers cited above, i.e., the 
redshift evolution of galaxies is not taken into account. Therefore, 
we will limit ourselves to comparison of the general picture of the 
S FR-dependence of the mass - metallicity relation. But first we 
review the results obtained by previous investigators. 

Ellison et al.l OOOS) considered a mass - metallicity relation 
for SDSS galaxies, binned in specific star for mation rate. Metal- 
licities were obtained using the calibration from lKewlev & Dopital 
( l2002h . They found that at high stellar masses (logM^ > 10) the 
mass - metallicity relation exhibits no dependence on the sS FR. 
At lower stellar masses, there is a tendency for galaxies with higher 
sS FR to have lower metallicities for a given stellar mass. The off- 
set in metallicity from highest to lowest sS FR bins is ~ 0.10-0.15 
dex in the stellar mass range 9 < logM, < 10. 

IMannucci et al. I (2010) studied the relation between stellar 
mass, metallicity and star formation rate for SDSS galaxies with 
redshifts 0.07 < z < 0.30. The oxygen abundances were derived 
using calibrations based either on the [N n]/i6584/HQ' or on R23 = 
([Oii]/13727 + [Oin]^4959,5007)/H8 parameter. When both cali- 
brations can be used, galaxy metallicity is defined as the average of 
these two values. Moreover, they selected only galaxies where the 
two values of metallicity differ no more than 0.25 dex. They found 
that galaxies with high stellar masses (logM, > 10.9) show no cor- 
relation between metallicity and star formation rate, while at low 
stellar masses more active galaxies also show lower metallicity. 

lYates et al. I ( 120121) studied the relation between stellar mass, 
metallicity and star formation rate for SDSS galaxies with redshifts 
0.005 < z < 0.25. Two metallicity-data sets were used. The oxy- 
gen abundances were derived following IMannucci et al. 1 ( 120100 in 
the first case and using the catalogued values in SDSS-DR7 in the 
second case. They found that the metallicity is dependent on star 
formation rate for a fixed mass, but that the trend is opposite for 
galaxies in the low and high stellar-mass ends. Low-mass galaxies 
with high 5 FR are more metal poor than quiescent low-mass galax- 
ies. Massive gal axies have lower metallic ity if their star formation 
rates are small. iLara-Lopez et al. I (l2013b have reached a similar 
conclusion. 

lAndrews & Martini I ( l2013n considered SDSS galaxies with 
redshifts 0.027 < z < 0.25. It has been assumed that galaxies at 
a given stellar mass, or simultaneously a given stellar mass and star 
formation rate, have the same metallicity. The SDSS galaxies were 
stacked in bins of (1) stellar mass and (2) both stellar mass and star 
formation rate. The auroral lines were measured in stacked spec- 
tra and the abundances were derived using the direct method. They 
found that the 5 Fi?-dependence appeared monotonic with stellar 
mass. 

In summary, the general conclusion in previous studies, is 
that the mass - metallicity relation is dep endent of (specific) 
star formation r ate at low stellar masses (lEllison et al.l 120081: 
'Mannucci et al. ' '20ltf; Yates et al.l I2OI2I : lAndrews & Martini I 
12013: Lara-Lopez et al. ii2013l) . Our results for low-mass galaxies 
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are in agreement with tiiese results. The conclusions by various au- 
thors regarding the 5 F/?-dependence of the mass - metallicity re- 
lation ayiigh_stellar masses ar e controversial, however . On the one 
hand jEUison et alj i200§) and lMannucci et al.l ^20 id) have found 
that at high stellar masses the mass - metallicity relation shows no 
dependence on (specific) star formation rate. Our res ults for mas- 
sive galaxies are in a greemen t with the conclusions by Ell ison et al 



2008) and Mannucci et al. 1 12010). On the other hand, Yates et al 



2012l) : lAndrews & Martini] ( l2013h : lLara-L6pez et"an ( l2013h have 
found that at high stellar masses the mass - metallicity relation 
depends on (specific) star formation rate. Our results for massive 
galaxies are in conflict with those results. 



4 LOW-MASS GALAXIES: SPIRALS VS IRREGULARS 

Thus, there is no coiTelation between oxygen abundance and sS FR 
in massive galaxies and the oxygen abundance correlates with the 
sS FR in low-mass galaxies. The difference in behaviour of massive 
and low-mass galaxies may be caused by the fact that our sample 
of low-mass galaxies is a mixture of spiral and irregular galaxies 
while massive galaxies can be expected to constitute subsamples 
of galaxies of a single morphological type (spiral galaxies). Below, 
we compare the properties of irregular and spiral galaxies of low 
masses. 



4.1 Spirals versus irregulars 

Let us firstly test whether the spiral and irregular galaxies of sim- 
ilar masses from our sample have different oxygen abundances. 
Morphological classifications of the SDSS galaxies have been per- 
formed within the framework of "Galaxy Zoo" projecjj. The data 
catalogue are available i n SDS S DR9 database. The project is de- 
scribed in Lintott et al. I ( 120081) and the data release is described in 
iLintott et al. I ( 1201 If) . The fraction of the votes in each of the six 
categories (elliptical, clockwise spiral, anticlockwise spiral, edge- 
on disc, don't know, merger) is given, along with "debiased" votes 
in elliptical and spiral categories and flags identifying systems as 
classified as spiral, elliptical or uncertain. We assume the votes of 
spiral galaxies is the sum of the fraction of the votes of clockwise 
spiral and counterclockwise spiral. Thus, each galaxy is described 
by five values fx which specify the probability that the galaxy be- 
longs to one of those classes. We assume that the galaxy is a spiral 
galaxy if the value of fsp > 2/3. Since the class "irregular galaxies" 
is not included, we select iiTegular galaxies by elimination, i.e., we 
assume that the galaxy is an irregular galaxy if each value fs,,, /e, 
/Edge, fuerg IS less than 1/3. 

If there is a difference between oxygen abundances in irregular 
and spiral galaxies of similar masses then one may expect the mass 

- metallicity relation for irregular galaxies to be different from that 
for spiral galaxies. The panel al) of Fig.|4]shows the mass - metal- 
licity diagram for spiral galaxies at the present epoch (for galaxies 
with redshifts less than 0.05). Light-gray (blue) points are individ- 
ual galaxies. Oxygen abundances in massive galaxies do not show 
significant correlation with galaxy mass. This flattening of the mass 

- metallicit y relation (a plateau at h i gh masses) has prev iously been 
observed bv lTremonti et al. l ( l2004 ; iThuan et an ( l2010l) . The lumi- 
nosity - metallicity relation f or nearby spiral galax ies also exhibits 
a plateau at high luminosities dPilvugin et al. 120070 . Since our goal 



is to compare the mass - metallicity relations for spiral and irregu- 
lar galaxies we should derive the relations for a galaxy mass range 
well populated with both spiral and iiTegular galaxies, logM^ < 
10.25. The best linear least-squares fit to the data for spiral galax- 
ies with masses 10.25 > logM^ > 9.0 is 

12 + log(0/H) = 6.20(±0.05) + 0.233(±0.006) x logM^, (14) 

shown in the panel al) of Fig. |4] as the dark-gray (red) solid line. 
The small formal errors in the coefficients of Ea.(ll4t and the rather 
small mean deviation (0.057 dex) are due to the relatively large 
number of data points (1494). 

The panel a2) of Fig. |4] shows the mass - metallicity diagram 
for irregular galaxies at the present epoch, z g 0.05. Light-gray 
(green) points are individual galaxies. The linear best fit to the data 
for irregular galaxies with masses 10.25 > logM^- > 8.5 



12 + log(0/H) = 5.68(±0.10) - 0.282(±0.01 1) x log Ms 



(15) 



is presented in the panel a2) of Fig.|4]by the dark (black) dashed 
line, while the dark-gray (red) solid line in the lower panel is the 
same as in the panel al). Comparison of these slopes, as well as 
comparison between Eq.(ll4t and Ea.( ll5t . shows that the mass - 
metallicity relation for iiTegular galaxies differs slightly from that 
for spiral galaxies. Note in passing, however, that the difference is 
rather small, less than the scatter in abundances in both irregular 
and spiral galaxies. 

Panel bl) of Fig. |4] shows the mass - sS FR diagram for spiral 
galaxies with redshifts less than 0.05. Light-gray (blue) points are 
individual galaxies. The best linear least-squares fit to the data for 
spiral galaxies with masses 10.25 > logM^- > 9.0 is 



9 + \og(sSFR) = 0.81(±0.28) - 0.17(+0.03) x logM^-, 



(16) 



The catalogs are available at|http://www.galaxyzoo.org] 



shown in panel hi) of Fig.|4]as the dark-gray (red) solid line. Panel 
b2) of Fig.|4]shows the mass ~ i'5 FR diagram for irregular galaxies. 
The linear best fit to the data for irregular galaxies with masses 
10.25 > logMi- > 8.5 

9 + \og(sSFR) = 4.51(±0.43) - 0.55(±0.05) x logMs, (17) 

is presented in panel b2) of Fig. |4] as a dark (black) dashed line, 
while the dark-gray (red) solid line in the lower panel is the same 
as in panel bl). Comparison of Eq.(fT6t and Eq.(|17K shows that the 
mass - sS FR relation for irregular galaxies is steeper than that for 
spiral galaxies. Spiral and irregular galaxies with masses logM^- ~ 

10 have similar i'5 FRs, while irregular galaxies with masses logM^- 
~ 9 have higher i'5 FRs than spiral galaxies of similar masses. 

Panel cl) of Fig. |4] shows the sS FR - metallicity diagram 
for spiral galaxies with redshifts less than 0.05. Light-gray (blue) 
points are individual galaxies. The best linear least-squares fit to 
the data for spiral galaxies with masses 10.25 > logM^ > 9.0 is 

12+log(0/H) = 8.44(±0.01)-0.058(±0.007)x(log(55Fi?)+9),(18) 

shown in panel cl) of Fig.|4]as a dark-gray (red) solid line. Panel 
c2) of Fig. |4] shows the i5 FR - metallicity diagram for irregular 
galaxies. The linear best fit to the data for irregular galaxies with 
masses 10.25 > logM^ > 8.5 

12-Hlog(0/H) = 8.23(±0.01)-0.160(±0.011)x(log(s5F/?)-^9),(19) 

is presented in panel c2) of Fig. |4] as a dark (black) dashed line, 
while the dark-gray (red) solid line in the lower panel is the same 
as in the panel cl). Comparison of Eq.(IT8t and Eq.(|19t, shows that 
the i5 FR - metallicity relation for irregular galaxies is steeper than 
that for spiral galaxies. 

The fact that our sample of low-mass galaxies is a mixture 
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Figure 4. Panel al). The mass - metallicity diagram for spiral galaxies. Light-gray (blue) points show individual objects, while the dark-gray (red) solid line 
shows the best fit to galaxies with the masses between \ogMs = 9.0 and logM^ = 10.25. Panel a2). The mass - metalhcity relation for irregular galaxies. 
Light-gray (green) points are individual objects, while the dark (black) dashed line is the best fit to galaxies with masses between logM^- = 8.5 and logMs 
= 10.25. The dark-gray (red) solid line is the same as in panel al). Panels bl) and h2) show the mass - sSFR diagrams for the same samples of spiral and 
irregular galaxies, respectively. Panels cl) and c2) show the sSFR - metallicity diagrams for spiral and irregular galaxies, respectively. (A colour version of 
this figure is available in the online version.) 



of spiral and irregular galaxies may explain the dependence of 
the metallicity - redshift relation on the sS FR observed for low- 
mass SDSS galaxies. The subsample of low-mass galaxies with star 
formation rates close to sS FR„,„y contains mainly irregular galax- 
ies while the subsample of low-mass galaxies with star formation 
rates around Q.lxsS FR,„a., contains both irregular and spiral galax- 
ies. Therefore these subsamples show different properties and the 
metallicity - redshift relations for low-mass galaxies with sS FR = 
sSFR^ax and sS FR = O.lxsS FR„„y are different. 



It is expected that in general the oxygen abundance in a galaxy 
correlates with three parameters, i.e., O/H = f(Ms , sS FR, z). Using 
the same data points as in the case of Ea.l ll4t . we have found for 



spiral galaxies 
12 + log(0/H) 



6.18 (±0.05) -1-0.240 (±0.006) logM^ 
0.023 (±0.005) (log sS FR + 9) (20) 

1.97 (±0.22) z. 



Again, using the same data points as in the case of Eq.ljTs), we have 
found for irregular galaxies 

12 -1- log(0/H) = 5.90(±0.11)-i-0.262(±0.001)logMs 

- 0.069(±0.009)(log.s5Fie-h9) (21) 

- 2.13(±0.56)z. 

Comparison of the coefficients in Ea.ll20t and Eq. Olt shows that 
the largest difference (by a factor of ~3) is found for the coefficients 
in the term with sS FR. 
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Figure 5. Upper panel. The mass - metallicity relation for irregular galax- 
ies. The gray (green) points show SDSS galaxies (same as in Fig.|4), while 
dark (black) points show the best fit to data for stellar masses between 
logA/i' = 8.5 and logM^- = 10.25. The dark (black) sohd line is the relation 
for nearby irregular galaxies from Berg et al. (2012) and the cr osses show 
the ma ss - metallicity relation for nearby irregulars derived by ILee et al. I 
12009) with constant mass-to-light ratios, while the plus signs show their re- 
lation for Mi- obtained with colour-dependent mass-to-light ratios. The dark 
(black) dashed line is the mass - metallicity relation for SDSS galaxies from 
IZahid et al. I j2012l) . Lower panel. The luminosity - metallicity relation for 
irregular galaxies. The gray (green) points are the SDSS galaxies, while the 
dark (black) points show the best fit to data for luminosities between Mg = 
-20 and Mg = -17. The dark (b lack) sohd hne is the relation for nearby irreg- 
ular galaxies from lBerg et al. (2012) and the crosses show the lumino sity - 
metallicity relation for nearby irregulars according to Lee et al. I I2006h . The 
dark (black) dashed l ine is the Mb - (O/H)?-^ relation for nearby irregulars 
from lPilvugin et al. I J2004h . (A colour version of this figure is available in 
the online version.) 



4.2 Comparison of the mass - metallicity and luminosity - 
metallicity relations for SDSS and nearby irregulars 

4.2.1 The mass — metallicity relations 

Fig.|5]shows a comparison between our mass - metallicity and lu- 
minosity - metallicity relations for SDSS galaxies and relations 
obtained for nearby irregular galaxies. The upper panel shows the 
mass - metallicity diagram, where the gray (green) points are irreg- 
ular galaxies from the SDSS, the same as in Fig.|4] The dark (black) 
points show the best fit t o the data in the mass range from logMs 
= 8.5 to logM.s - 10.25. iLee et al. I ( |2006) have constructed mass 



- metallicity and luminosity - metallicity relations for 25 nearby 
irregular galaxies. The oxygen abundances in 21 of these galaxies 
were derived through the T^ method. They estimated the masses of 
galaxies from their luminosities in two ways. In the first case they 
adopted a constant stellar-mass-to-NIR-light ratio and obtained the 
following mass - metallicity relation, 



12 + log(0/H) = 5.26(±0.27) + 0.332(±0.033) x log Ms 



(22) 



This relation is shown in the upper panel of Fig. |5]by dark (black) 
crosses. In the second case they adopted a colour-dependent stellar- 
mass-to-NIR-light ratio, which yielded the following mass - metal- 
licity relation. 



12 + log(0/H) = 5.65(±0.23) + 0.298(±0.030) x logM^ 



(23) 



This relation is shown in the upper panel of Fig. |5]by dark (black) 
plus signs. Furthermore, iBerg et al. I (120121) measured the temper- 
ature sensitive [Om]/i4363 line for 31 low luminosity galaxies in 
the Spitzer Local Volume Legacy survey, and thus determined oxy- 
gen abundances using the direct method. They have analysed a 
"Combined Select" sample composed of 38 objects (taken from 
their sample and as well as the literature) with direct oxygen abun- 
dances and reliable distance determinations. Stellar masses where 
estimated from the 4.5 yum luminosities and K ~ [4.5] and B ~ K 
colours. Using these data, they derived a mass - metallicity rela- 
tion, 



12 + log(0/H) = 5.61(±0.24) + 0.29(±0.03) x XogMg 



(24) 



This relation is shown as the dark (black) solid line in the upper 
panel of Fig. [5] 

The upper panel of Fig. |5] shows that our mass - metallic- 
ity relation, (0/H)co;^ - Ms , for SDSS galaxies are consistent with 
(0/H)7-, - Mg re latio ns for nearby irre gular galaxies obtained by 
iLee et al. I (|2006|) andlBerg et al. I J20TI) . 

Zahid et al. I ( 120121) have estimated oxygen abunda nces in the 



SDSS galaxies using the empirical N2 calibration after I Yin et al. I 
( 120071) and have constructed a mass - metallicity relation, 

12 -I- log(0/H) = 5.585(±0.003) + 0.47(±0.01) x (log Ms - 9)(25) 

which is shown in the upper panel of Fig.|5]by dark (black) dashed 
line. This mass - metallicity relation for SDSS galaxies differs sig- 
nificantly from the mass - metallicity relation for SDSS galax- 
ies considered here an d from the t he ma ss - metallicity r e lation s 
for nearby galaxies by iLee et al. I ( 20061) and Berg et al. I (120121) . 
This may appear suiprising, since 



Zahid et al. 



(2012) have esti- 



mated oxygen abundance s using the purely empirical N2 calibra- 
tion. However, lYin et al. I ( 120071) have noted that the N2 index is 
only a useful metallicity indicator for galaxies with 12-1- log(0/H) 
< 8.5 as they had no calibrating data points with 12 + log(0/H) > 
8.5. Hence, the high oxygen abunda nces (up to 12 + log(0/H) ~ 
9.0) obtained by IZahid et al. I ( 120121) with this calibration may be 
disputed. 



4.2.2 The luminosity - metallicity relations 

The lower panel of Fig. |5] shows the luminosity - metallicity dia- 
gram for irregular galaxies. The gray (green) points show 0/H - 
Mg diagram for irregular galaxies from our SDSS sample, where 
Mg is the absolute magnitude of a galaxy in the SDSS photometric 
band g. The dark (black) points show the best fit to those data. 



12 -I- log(0/H) = 6.41(±0.16) - 0.103(+0.009) x Mg 



(26) 
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obtained in tiie luminosity range from Mg = -20 to Mg = -17. The 
B - g colour index is ~ 0. 1 mag only for the star-forming galaxies 
and t he 0/H - Mg and O/H - Mb relations are t herefore compa- 



rable jPapaderos et al. ll2008HGuseva et al. |[2009l) . The O/H - Mg 
relation obtained bv lSerg et al. 1(120121) for their "Combined Select" 
sample of nearby irregular galaxies, 



12 + log(0/H) = 6.27(±0.21) - 0.11(±0.01) x Mg, 



(27) 



is shown by the dark (black) solid line. It s hould be noted that thi s 

relation seems to be valid up to Mg 8 (Skillman et al. 20131). 

The relation for nearby irregular galaxies by Le e et al. (2006 ), 

12 + log(0/H) = 5.94(+0.27) - 0.128(+0.017) x Mb, (28) 

is shown by dark (black) crosses, and that bv lPilvugin et al. I ( |2004) . 

12 + log(0/H) = 5.80(+0.17) - 0.139(+0.011) x Mb, (29) 



is represented by the dark (black) dashed line. 

Our luminosity - metallicity relation, (O/H)^-,^ - Mg 



for 



SDSS galaxies is in satisfactory agreement with (0/H)7, - Mb rela- 
tions for nearby irregular galaxies obtained in several oth er studies 
teerg et al. ll201^ : lLee et 571120061 : IPiTvugin et al. II2004I) . It should 
also be noted that the galaxy metallicity correlates more tightly 
with its stellar mass than with its luminosity. 



5 DISCUSSION 

5.1 Does the aperture effect contribute to the redshift 
dependence of the oxygen abundances? 

The right panels in Fig.|3]show there is a redshift dependence of the 
oxygen abundance in galaxies of different masses in the sense that 
the oxygen abundances in distant galaxies are lower than that in 
nearby galaxies. Chemical evolution of galaxies is usually consid- 
ered the only reason for the redshift variation of the oxygen abun- 
dance in galaxies. In the case of the SDSS galaxies there seems to 
be another effect that could make a contribution to this redshift de- 
pendence. The SDSS galaxy spectra span a relatively wide range 
of redshifts. There is thus an aperture-redshift effect present in the 
SDSS data (all spectra are obtained with 3-arcsec-diameter fibers). 
At a low redshift the projected aperture diameter is smaller than that 
at high redshift. This means that, at low redshifts, SDSS spectra can 
be seen as global spectra of the central parts of galaxies (if the fibers 
are centered on centers of galaxies), i.e., composite nebulae includ- 
ing Hii regions near the center. At large redshifts, however, SDSS 
spectra are closer to global spectra of whole galaxies, i.e., spectra 
of composite nebulae including H n regions near the center as well 
as in the periphery. Since there is u sually a radial abundance gradi- 
ent in the discs of spiral galaxies I Vila-Costas & Edmunds 1119921 : 



IZaritskv et al. lll994l : |PilYugin et al. l2004l among others), the oxy- 
gen abundance in the composite nebulae in distant galaxies may 
be lower than that in nearby galaxies if the fibers are centered on 
centers of galaxies. 

We estimate the aperture-correction factor A (see Eq|12t to 
account for the fact that only a limited amount of emission from a 
galaxy is detected through the 3" diameter fiber The quantity 1/A 
is the galaxy light fraction within the fiber The left panels of Fig.|6] 
show the redshift dependence of the 1/A for galaxies of different 
masses. The right panels of Fig. |6] show the oxygen abundance in 
a galaxy as a function of 1/A for subsamples of galaxies of differ- 
ent masses. Light gray (blue) points a show individual objects with 
redshifts in the full considered range, 0.02 < z < 0.33 and the dark 



Table 1. Coefficients in the "oxygen abundance - redshift" relations, 12 + 
log(0/H) = a + hxz, for subsamples of galaxies of different stellai' masses 
and different fraction of the light within fiber 1/A. (A colour version of this 
figure is available in the online version.) 



mass range 
logMi- 



range of 1/A 



10.25 - 10.75 0.45 > 1/A > 0.35 8.591±0.004 -0.456±0.037 

1 > 1/A > 8.596±0.001 -0.488±0.012 

10.75-11.25 0.45 > 1/A > 0.35 8.615±0.005 -0.326±0.034 

1 > 1/A > 8.615±0.002 -0.325±0.012 

11.25-11.75 0.45 > 1/A > 0.35 8.623±0.012 -0.250±0.057 

1 > 1/A > 8.617±0.004 -0.222±0.022 



(black) line is the best fit to the data. Dark gray (red) points show 
individual objects with redshifts in the range 0.08 - 0.1. There is 
obviously a correlation between the oxygen abundance in a galaxy 
and the galaxy light fraction within the fiber. The correlations for 
the full considered redshift range and that for a restricted redshift 
range are similar. The galaxy light fraction within the fiber also cor- 
relates with the redshift. Therefore, the diagrams in Fig. |6] cannot 
tell us whether the correlation between the oxygen abundance in a 
galaxy and the galaxy light fraction within the fiber is caused by 
redshift evolution (emichment in heavy elements) of galaxies or by 
the aperture effect or a combination of both. 

To solve this dilemma we compare the redshift variation of the 
oxygen abundance in subsample of galaxies for a given value of the 
aperture correction factor and that for total sample of galaxies of a 
given mass. Gray (blue) points in the upper panel of Fig.|7]show the 
redshift dependence of the oxygen abundance with galaxy masses 
logMs ranging from 10.25 to 10.75, and 1/A in the range 0.45 > 
1/A > 0.35. The dark (black) crosses is the linear best fit to those 
data 



12 + log(0/H) = a + bxz. 



(30) 



Values for the coefficients a and b of Eq.OQIl are listed in Table [T] 
The small statistical uncertainties in the coefficients a and b is due 
to the large numbers of galaxies in the subsamples. The dark (black) 
line is the linear best fit to all galaxies within the mass range con- 
sidered, which essentially coincide with the crosses representing 
the subsample. 

Similar diagrams for the other mass intervals are presented in 
the middle and lower panels of Fig. |7] Again, the relation O/H - z 
for the subsample and the total sample of galaxies within the given 
mass ranges coincide. This means that the aperture effect does not 
make any significant contribution to the redshift dependence of the 
oxygen abundance in galaxies. The trend in oxygen abundance with 
redshift must be caused by the chemical evolution of galaxies, i.e. 
by the increasing level of astration at low redshifts. 

The fact that the aperture effect does not play any significant 
role in the redshift dependence of the oxygen abundance may be 
explained if the fiber positions are distributed randomly over the 
images of galaxies. If this is the case, even in nearby galaxies, the 
fibers may cover both the central and the periphery parts of the 
galaxies. In such case, the average oxygen abundances in a subsam- 
ple of SDSS galaxies is more similar to the characteristic oxygen 
abundance in a galaxy rather than to its central oxygen abundance. 
It should be noted that many observations with the fibers centered 
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Figure 6. Left panels. Redshift dependence of the galaxy light fraction within the fiber (1/A) for galaxies of different masses. Right panels. 
abundance in a galaxy as a function of 1/A. Light gray (blue) points are individual objects with redshifts from the full considered range, 0.02 < z 
(black) hne is the best fit to those data. Dark gray (red) points are individual objects with redshifts from the range 0.08 - 0.1. (A colour version 
is available in the online version.) 



The oxygen 
< 0.33. Dark 
of this figure 



on the galaxy nuclei are excluded from our analysis since the spec- 
tra are often classified as AGNs and not star- forming regions. 

5.2 Time-averaged specific star formation rate in massive 
galaxies 

The trend in oxygen abundance with redshift is caused by the in- 
creasing level of astration (or by the decreasing of the gas-mass 
fraction jx = MgaJiMs+Mg^,)) with decreasing of redshift. This 
provides a possibility to estimate the value of the star formation 
rate in galaxies of different masses averaged over a fixed time in- 
terval. We will consider here only massive galaxies, which can be 
expected to constitute subsamples of galaxies of a single moipho- 
logical type (spiral galaxies). 

The prediction of the simple (closed-box) model of galactic 
chemical evolution can be well approximated by a linear expression 
for /J values between ~ 0.7 and ~ 0.1. In a real situation, the oxygen 
abundance is also afi' ected by the m ass exchange between a galaxy 
and its environment ( |Pagellll997D . It is commonly accepted that 
gas infall plays an important role in the chemical evolution of disks 
of spiral galaxies. Therefore, the application of the simple model to 
large spiral galaxies may appear unjustified. It is expected that the 
rate of gas infall on to the disk decreases expo nentially with time 
|Calura et al. Il2009l e.g.). It has been shown iPilvugin & Ferrini I 
Il998j ) that the present-day location of a system in the yu - O/H di- 
agram is governed mainly by its evolution in the recent past, but 



is weakly dependent of its evolution on long timescales. Therefore, 
the fact that the present-day location of spiral galaxies is near that 
predicted by the simple model is not in conflict with infall of gas 
onto the disk taking place over a long time (the latter is assumed 
necessary to satisfy the observed abundance distribution function 
and the age - metallicity relation in the solar neighbourhood) since 
these observational data reflect the evolution of the system in a dis- 
tant past. A decrease of fj by O.I results in a i ncrease of oxygen 
abundance by ~ 0.13 dex jPilvugin et al. Il2006h . which lead us to 
the approximation 



A(log(0/H))a;-1.3xA;U. 



(31) 



The variation in oxygen abundance with redshift is given by the 
cooefficient b in Eq. OOb . the value of -0.1b corresponds to a 
change in oxygen abundance A(log(0/H)) over a redshift interval 
Az = O.I. At low redshifts, a redshift interval Az = 0. 1 corresponds 
to a time interval of about 1 Gyr. Taking Eq. (l3It into account, the 
change of gas-mass fraction A/i over the time interval of 1 Gyr cor- 
responds to Q.Qllb. The gas mass fraction is expressed in terms of 
the total mass of a galaxy M,o,. The specific star formation rate is 
expressed in terms of the stellar mass of a galaxy Ms . The oxygen 
abundance 12 + log(0/H) = 8.6 (which is a typical value for mas- 
sive gal axies, Fig.[7]> correspo nds to the value of gas mass fraction 
H ~ 0.2 jPiIvuginetani2006h . This gives Mg ^ 0.8x/W„„. Then the 
change of gas mass per year expressed in terms of the stellar mass 
of a galaxy (i.e. the average sS FR) is given by the value lO^'^xfo. 
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Figure 7. Redshift dependence of the oxygen abundances in galaxies of 
different masses and an aperture correction factor 0.45 > 1/A > 0.35. Gray 
(blue) points are individual objects, while the dark (black) crosses show the 
best fit to the data. The dark (black) fine is the best fit to all galaxies in the 
given mass range. (A colour version of this figure is available in the online 
version.) 



Using the values of b from the Table [T] we estimate the time- 
averaged sS FRs for massive galaxies. We find 9+log(sS FRgy^,) ~ 
-1.31 for galaxies of masses logM^- = 10.5 +/-0.2, 9+log(sS FRaver) 
~ -1.49 for galaxies of masses logMi- = 11.0 +/-0.2, and 
9+log(i'5 FRaver) — 1-64 for galaxies of masses logMs = 11.5 +/- 
0.2. Comparison between the 55 FR,„„y (Fig.[3]( and sS FRi,y„ shows 
that the sS FR,„ax exceeds the sS FRaver by a factor of 6-7. Since the 
sS FRmax we obtain is not an absolute maximum star formation rate 
but the mean value for only 5-10% of galaxies with highest sS FRs 
(see Fig.[3]l, the absolute maximum star formation rate can exceed 
the time-averaged star formation rate in galaxies of a given mass 
by about an order of magnitude. 



counterpart. In total, we used a sample of 233 reference Hit re- 
gions for abundance determinations. The conditions |log(0/H)c„j - 
log(0/H)coJ < 0.05 and |log(N/H)c;,, -log(N/H)coJ < 0.10 were 
used to select the SDSS objects with the most reliable abundances. 
Only the (O/H)^^^ abundances were used in our further analyses. 

We analysed the redshift dependence of the maximum specific 
star formation rate sS FRi„ay in SDSS galaxies of different masses. 
Six stellar mass intervals centered on the values logM., = 9.0, 9.5, 
10.0, 10.5, 11.0, 11.5 and width interval widths +/- 0.2 dex were 
considered. We find that the sS FR^ax decreases with increasing 
stellar mass and decreasing redshift. For each mass interval we 
considered the metallicity - redshift relations for subsamples of 
galaxies with sS FR ~ sSFR„,a-, and with sS FR ~ O.lxsS FR„,„^. 
These relations coincide for massive (logM^- > 10.5) galaxies, i.e., 
the metallicity ~ redshift relations for massive galaxies are inde- 
pendent of the i'5 FR. On the other hand, these relations differ for 
low-mass galaxies as the galaxies with sS FR ~ sSFR„,ax have, 
on average, lower oxygen abundances than galaxies with sS FR ~ 
OAxsS FR,„a,. 

We also find evidence in favour of irregular galaxies having 
higher (on average) sS FRs and lower oxygen abundances than 
spiral galaxies of similar masses and that the mass - metallicity 
relation for spiral galaxies differs slightly from that for irregular 
galaxies. The mass - metallicity and luminosity - metallicity rela- 
tions obtained for irregular SDSS galaxies agree with the same for 
nearby irregular galaxies with direct abundance determinations. 

The fact that any sample of low-mass galaxies is likely a mix- 
ture of spiral and irregular galaxies can explain the dependence of 
the metallicity - redshift relation on the i'5 FR. 

The redshift variation of the oxygen abundances in galaxies of 
a given mass allows us to estimate the decrease of gas-mass frac- 
tions, i.e. this provides a possibility to estimate the star formation 
rate in galaxies of different masses averaged over a given time inter- 
val. We estimated the time-average specific star formation rates for 
massive galaxies. We find that the absolute maximum star forma- 
tion rate can exceed the time-averaged star formation rate in mas- 
sive galaxies by about an order of magnitude. 

Finally, we compared the relations O/H - z for subsample of 
galaxies within a given range of galaxy light fractions within the 
spectroscopic fiber and for total sample of galaxies of a given mass. 
These subsamples show very similar relations, which suggests the 
aperture effect does not make a significant contribution to the red- 
shift dependence of the oxygen abundances in galaxies. This may 
be explained by fiber positions being distributed randomly over the 
images of the galaxies, so in galaxies the fibers may cover both cen- 
tral and periphery parts of different galaxies. In such case, the aver- 
age oxygen abundances in a subsample of SDSS galaxies is more 
similar to a characteristic oxygen abundance in a galaxy rather then 
to its central value. 



6 SUMMARY AND CONCLUSIONS 

In the present study, we have used a sample of 98.986 star-forming 
SDSS galaxies to examine the dependence of the metallicity - red- 
shift relations on the star formation rate for galaxies of different 
stellar masses. We determined two values of oxygen and nitro- 
gen abundances for each object using a modified version o f the 
recent counterpart method (C method) jPilvugin et al. l2012a) . The 
(0/H)coN and (N/H)coN abundances were derived using the logi?3, 
P and logiNiIRi) values to find a counteipart for the studied SDSS 
galaxy, and, similarly, the (0/H)cns and (N/H)cns abundances were 
derived using the logi?3, \0gN2 and log(A'2/52) values to find a 
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